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Abstract 

We follow a low-energy effective theory approach to identify the general class of theories that 
describes a real vector field (of unconstrained norm) coupled to gravity. The resulting set may be 
regarded as a generalization of the conventional vector- tensor theories, and as a high-momentum 
completion of aether models. We study the conditions that a viable cosmology, Newtonian limit 
and absence of classical and quantum instabilities impose on the parameters of our class of models, 
and compare these constraints with those derived in previously studied and related cases. The 
most stringent conditions arise from the quantum stability of the theory, which allows dynamical 
cosmological solutions only for a non-Maxwellian kinetic term. The gravitational constant in the 
Newtonian limit turns to be scale dependent, suggesting connections to dark matter and degravi- 
tation. This class of theories has a very rich gravitational phenomenology, and offers an ample but 
simple testing ground to study modifications of gravity and their cosmological implications. 



I. INTRODUCTION 



The symmetries of a spacetime severely constrain the transformation properties of its 
matter content. In Minkowski spacetime for instance Lorentz symmetry allows scalar fields 
to have a non-vanishing expectation value, forbidding non-scalar fields from being non- 
zero. It is perhaps because of these constraints that most cosmologists have focused on 
cosmological solutions with scalar fields. 

Nevertheless, in an expanding universe no symmetry principle prevents the existence of 
non-vanishing tensor fields, as long as they are constant in space and invariant under spatial 
rotations. Since only bosonic fields can acquire large occupation numbers, and if we restrict 
ourselves to fields of spin less than two, this leaves vector fields (or one-forms) as essentially 
the only additional type of field that may be relevant cosmologically. If the vector field only 
couples to gravity, the presence of a preferred Lorentz frame will remain unobservable in the 
matter sector, while, of course, Lorentz symmetry is broken anyway by the expansion of the 
universe and is not of particular relevance on large scales. 

Massless vector-tensor theories were considered long time ago to explore different modi- 
fications of gravity, and as alternatives to the more prevalent scalar-tensor theories at that 
time [H E]. Indeed, as we shall see, vector-tensor theories have a much richer gravita- 
tional phenomenology than scalar-tensor theories, and thus offer a wider framework to 
study departures from general relativity. Theories with vectors have also experienced a 
modest revival in connection with inflation and the accelerated expansion of the universe 
[SliaElinilTllHlinilinillllliailHllTlllISllinillT^ studies of the spontaneous breaking of 
Lorentz invariance in the so-called Einstein-aether models of Jacobson and Mattingly [18\ 
(see fT9j for earlier discussions, and [20] for a status report.) 

In this article we study a class of theories that somewhat differs from previous models. 
We explore the most general vector-tensor theory expected to describe our universe at low 
energies, large distances and long timescales. We thus follow a low energy effective action 
approach and consider all generally covariant terms with the lowest number of vector field 
derivatives. Since we deal with massive vectors, a consistent quantum description at low 
energies does not require gauge invariance, so many different kinetic terms for the vector 
field are possible. Conversely, if we do not assume gauge invariance, no symmetry principle 
prevents the appearance of mass terms for the vector, or vector self-interactions for that 
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matter. Our effective action thus contains all generally covariant terms with two vector field 
derivatives and a general quartic potential^ V that general covariance forces to depend on 
the squared norm of the vector, A^A'^. 

Up to a constant term, any scalar quartic potential can be written in the form 

V{A^A^) = X{A^A^ + My, (1) 

where M is a mass scale and A a dimensionless parameter. In fact, around an extremum, 
any potential can be locally approximated by ([T]). Such a dependence opens up a quite 
interesting possibility, because the vacuum of the theory has a non-vanishing vector field 
A^A^^ = — M^, and thus in fiat space Lorentz symmetry is spontaneously broken.^ In fact, 
in some modern approaches to quantum gravity Lorentz invariance is expected to be broken 
at high energies [22] • However, if broken, it should be done in a very subtle way. If 
the breaking relies in the existence of a fundamental preferred frame at the Planck-scale, 
some unnatural fine tuning issues appear p3j, or additional ill-behaved modes appear in 
the spectrum [251 ES]. Moreover, a fixed external structure goes drastically against our 
currently accepted principles (i.e. general covariance). This suggests that, probably, the 
more reasonable way in which Lorentz invariance could be broken is by the existence of a 
vacuum expectation value (vev) for some non-scalar field, which, as in the ghost condensate 
[26] . may be the gradient of a scalar. 

A widely studied class of models in which Lorentz invariance is spontaneously broken 
involves the Einstein-aether models introduced by Jacobson and Mattingly in [18]. One 
of the most distinctive properties of these models is that the norm of the vector field is 
constrained to have a constant value. This constraint may appear artificial, but it could be 
naively justified as the limit in which the coupling constant A in the potential of equation 
([T]) tends to infinity. Hence, aether models may be regarded as a particular case of the class 
of theories we consider here. Because we drop the fixed- norm constraint we call this class 
of theories "unleashed aether" models. Similar kind of actions have been considered for 
instance in [23 [2H1 1211 EQ] under the name of "bumblebee models." In general, the latter 
are models in which a vector field acquires a vev due to the presence of a potential. In 

^ As we shall see, the "potential" V actually contains kinetic terms for some of the vector degrees of freedom. 
^ We consider here positive values of M^. In that case, is timclikc in the broken phase. If is negative, 
the vector vev is spacelike. 
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that sense, the unleashed aether is a bumblebee model. The real difference between the 
unleashed aether and the most commonly studied bumblebee models lies in the parameters 
in our action that are assumed to be zero in the bumblebee case. 

The main goal of this article is to study the conditions that self-consistency and phe- 
nomenological viability imposes on the parameters of our class of theories, and to determine 
to what extent these differ from massless vector fields and aether models. In Section [III we 



present and discuss general properties of our class of models. In Section III we study the 



properties of cosmological solutions at early and late times, and derive basic constraints on 



the value of the scale M in the potential. In Section IV we study the Newtonian limit of 
these theories, and explore their relation with massless theories and aether models. Finally, 
in Section |V| we derive the conditions that classical and quantum- mechanical stability places 
on the parameter space of unleashed aether models. 

II. GENERALITIES 

We want to identify the most general action that describes the dynamics of a real vector 
field coupled to gravity. Without additional restrictions this action will contain all 
possible terms compatible with the symmetries of the theory, which we assume to consist 
just of general covariance and a Z2 symmetry —>■ —A'^. If we are interested in low 
energies however, we expect only those terms with the least number of derivatives to be 
relevant in the limit of long distances and times. We can thus concentrate on the most 
general generally covariant action with at most two derivatives acting on the metric and 
the vector. In principle, such an action would also contain arbitrary powers of the vector 
field, but since a vector has dimensions of mass, terms with too many powers of A^ must 
be suppressed by corresponding powers of an energy scale A. If the background value of the 
field is smaller than A, the dominant terms are hence given by those with no suppression, 
that is, by the operators of dimension four. Therefore, we just have to consider the most 
general action with two derivatives acting on the metric and vector field, but with no vector 
field operators of dimension higher than four: 

R Pi 



S = d X \f—g 



\f^,F^- - fOy^A^f + (ii^R^^A^A- + 



(2) 
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where the j3i are dimensionless coefficients, F^i, = d^Ay—duA^, and is a quartic potential of 
the form ([T]), which we plot in Figure [l] This is the most general action with two derivatives 
acting on A^ because a term proportional to V^Aj.V'^A^ can be eliminated integrating 
by parts and using the identity [V^, Vjy]A^ = R^J^^jA'^ . The only dimension four operator 
with one vector derivative, A^A'^'V uA'^ , is excluded by the discrete symmetry A^ —A^. 
Observe that we do not assume gauge invariance, which would only allow the Maxwell term 
in the vector sector. Although gauge invariance seems to be necessary for a consistent 
description of massless vector field interactions [31], it is not required for massive vector 
fields, at least at sufficiently low energies. The stability of the theory and phenomenological 
considerations then set constraints on the values of the so far arbitrary coefficients jSi and 
the derivatives of the potential V. When we explore these constraints in Section |V| we shall 
see that the term proportional to the Maxwell Lagrangian describes a vector field (under 
Lorentz-transformations), whereas the term proportional to the divergence squared of the 
field actually introduces an additional scalar (under Lorentz transformations). In that sense, 
it is debatable whether this is actually just a vector field Lagrangian. 

Of course, we should also add to the action ^ all the matter terms of the standard model 
of particle physics. A general vector field action would contain couplings of the vector field 
A'^ to this sector, as in bumblebee models. Such "Lorentz-violating" interactions have been 
thoroughly considered and constrained in the Standard Model Extension program initiated 
by Kostelecky (see e.g. E3])- In this article we assume that such couplings are negligible. 
Since interactions between the vector field and the matter sector would be mediated by 
gravitation, it is consistent to assume that they are small. For a review of modern tests of 
Lorentz invariance, see |34j . 

The validity of the effective theory proposed here requires that the background value 
of the vector field components satisfy \A'^\ ^ A. This is the assumption that allows us to 
neglect terms in the Lagrangian such as, for instance, 

C D ^A'^A'^V.Af'VyA,. (3) 

Because for a potential like ([!]) we expect the field to be in the vicinity of A^A^^ = — M^, this 
implies that <^ A^. This raises a mild fine-tuning reminiscent of the hierarchy problem — 
why is the Higgs vev much smaller than the Planck scale. Up to an overall constant, a generic 



5 



FIG. 1: A plot of the potential ([l]) for > and A > 0. Because of the indefinite signature 
of spacetime, the potential does not have the shape of a Mexican hat. Instead, around the two 
local minima along the Aq direction, the potential looks like a horse saddle and thus has negative 
Gaussian curvature. If < the roles of ^-nd Ai are reversed. 

potential would be of the form 

V{A,A^') = 2\^A^ A,A^ + ^A^A^' + ^{A.A^' + ■■■, (4) 

with dimensionless coefficients Aj of order one. Comparing with equation ([T]) we find that 
A = A2 and = (Ai/A2)A^ Thus, < requires Ai < A2. This fine-tuning does not 
have to be severe though. It suffices for instance that A1/A2 = M^/A^ < 10~^ in order for 
our expression to be relatively accurate. On a more practical side, the restriction M < A 
is almost inescapable: If M were of the order of A, we would have to include all powers 
of A'^ in our effective theory. Incidentally, a very small value of A can be easily justified if 
P13 = /?4 = 0. In that case, as we discuss in Section |V} the theory has additional symmetries 
in the limit A 0. For a review about effective field theories see [S3] . 

The action ^ is to some extent a generalization of previously considered theories. Setting 
V = one recovers the class of vector-tensor theories studied in [36], while replacing the 
potential ([T]) by a fixed-norm constraint leads to a subclass of the Einstein-aether models 
considered in [18], and the vector field theories studied by Carroll and Lim ^37].^ Although 

^ Our parameters are related to those of [37] by /3i3 = Pi + P3, (3 = (3i + + P3, and by those of [36] by 
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a fixed normed constraint seems to be difficult to justify from a particle physics perspective, 
general potentials V{A^A'^) have been argued to appear in theories in which gauge invariance 
is spontaneously broken [21 |M1 ES], and a fixed- norm constraint can be naively recovered 
from a potential of the form ([T]), at least classically, by taking the limit A — » oo (we shall 
make this connection more explicit below.) 

We are ready now to explore some of the physical consequences of a vector field described 
by the action ([2]). Varying with respect to the field leads to the equations of motion 

/9iV,F% + pV,{V,An + ^ii:,R,,A^ + (3^RA^ - ^A, = 0. (5) 

Because the covariant derivative and the metric commute, the equations we obtain by 
varying with respect to A'^ are equivalent to the ones we get varying with respect to 
Afj_. In fact, because the transition from a vector to a form involves just a field redefi- 
nition, A'^ ^ A^ = g^yA'^ ^ assuming that the action ([2]) is a functional of a vector A^ or 
a form A^ leads to physically equivalent results both classically and at the quantum level, 
even though, curiously, the corresponding expressions for the energy-momentum tensors 
T^i, = —{2/^/—g){6SA/Sg'^'') have different forms: If the action is taken to be a functional 
of a form A^, the energy-momentum tensor is 

T';r = 2/3,F,,F,^ + 2/3 [(A''V,(V.A'^) + (V.A^^) g^, - 2A(^V,)(V,A0] + 
+ As [2VpV(^(A,)A'') - V,VM'An9^.u - n{A^A,) - AA^Rp^^A^^] + 
- 2/?4 [R^^uA.A^ + RA^A, - S/^VMpA") + U{ApA'')g^,] + 
+ 2^A^A, + Cg^,, (6) 

while if the action is a functional of a vector A^, the energy momentum tensor turns out to 
be 

+ApiRA^,A^ - A^A^A^. (7) 

The additional terms in the energy-momentum tensor of a vector are proportional to the 
equation of motion (|5]), and thus vanish identically. In the following, we work with the 
energy-momentum tensor of a form, equation rt6]). Note that dA'^ stands for d{Afj_A^). 



/3i = 2e — r, /3 = — r, P13 = i] — t and Pi = u. 



III. COSMOLOGICAL SOLUTIONS 



On large scales, our universe seems to be well described by a (spatially flat) Friedmann- 
Robertson- Walker (FRW) spacetime, 

ds^ = -dt^ + a?{t)dx'^. (8) 

Its symmetries, homogeneity and isotropy, already constrain the possible configurations of 
a non-vanishing cosmic vector field quite strongly. Homogeneity (invariance under spatial 
translations) implies that a cosmic vector field must be spatially constant, while isotropy 
(invariance under spatial rotations) implies that its spatial components must vanish.^ Hence, 
a cosmic vector field must be of the form 

A'^ = (n(t), 0,0,0). (9) 

Because of homogeneity and isotropy, the spatial components of the vector field ([9]) satisfy 
the equations of motion ([s]), while the /i = component of equation ([s]) dictates that u has 
to obey 

(3{u + SHii) + \v' - 3(/3i3 + 4f34)H^ + 3(/5 - As - 2/34)i^] u = 0, (10) 

where a dot denotes a derivative with respect to cosmic time t, a prime a derivative with 
respect to A^A^, and H = d/a is Hubble's constant. It follows automatically from equation 



(10) that if /3 = 0, M is non-dynamical. Indeed, in that case the "kinetic" term of the vector 
is the one of electromagnetism, which does not contain any time derivatives of Aq = —u. 
This is why there are no terms proportional to (3i in the equation of motion for u. In aether 
models, the fixed norm constraint simply demands that the value of u be constant. 

The energy-momentum tensor ^ also explains why isotropy requires that the vector field 
have zero spatial components. Otherwise, say, the term proportional to dV/dA^ would yield 
an energy-momentum tensor which would not be invariant under spatial rotations, and the 
FRW metric ([s]) would not be a solution of Einstein's equations. Certainly, if the vector 
field is subdominant, it is not strictly necessary that A'^ satisfy equation (|9]) in order for 
equation ([s]) to be an approximate solution of the Einstein equations. Ultimately, equation 
(|9| should be justified dynamically, though we shall not explore this issue here. If the vector 

^ Except for some especial cases in which non-homogeneous field configurations are pure gauge so do not 
contribute to the energy- momentum tensor. 
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does point in the time direction, the energy density p = — T°o and pressure pSj = T^j are 
given, from ([6]), by 



p = p{2uu - ii^) + 6(^3 + 2pi)Huu + 2V' + 6{p - P13 - 2pi)H - 9{p + 2pi)H 



p = 2{(3- f3i3 - 2/34) (mm + u') - f3u' + (3/3 - 2l3u - 213^) AHuii + {2H + 3H^)u'^ 



(11) 
-V. 

(12) 



All the non-diagonal components of the energy- momentum tensor vanish, as demanded by 
isotropy. Because the Maxwell term does not contain u{t), the energy density and pressure 
do not have factors proportional to Pi. 

It will be useful to work with a simplified form of the energy density. The equation of 



motion (10) can be used to eliminate the term proportional to u in p, leading to the simpler 



expression 

p = _ _ _ 2p^)Huu - 3(3/3 - 2/3i3 - 2P4)HV + V. (13) 



Using the equation of motion again, it can be easily shown that the energy density (13) 
satisfies the conservation equation p + 3H{p + p) =0, which is just an expression of the 
conservation law V^T'^j, = 0. Note that p may become negative for certain values of the 
parameters. In principle one should not be scared by the possible appearance of negative 
energies, as long as the Hamiltonian has a local minimum. We explore the classical and 
quantum stability of this class of theories in Section |V} The executive summary is: Unleashed 
aether theories are stable only if /3 = or /3i = 0. Since cosmological solutions with P = 
have no dynamics, we shall thus set Pi = 0.^ In any case, this does not have any effect 
on our equations, because, as mentioned above. Pi does not appear in the field equation 
of motion nor the energy-momentum tensor. For Pi ^ ^ P the theory is unstable, at 
least in the ultraviolet. However, by an appropriate choice of the parameters we can get the 
dangerous modes out of the spectrum, recovering a sensible low energy effective field theory: 
the aether model. 



^ It is important to point out here that we are only considering classical (coherent) excitations for the vector 
field. Theories with (3 — but (3i ^ could also have cosmological relevance (as it happens for instance 
with the electromagnetic field) if they do not behave classically. (For instance, in the case of the Cosmic 
Microwave Background, what we have is a thermal distribution of photons). 
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A. Solutions at early times 



We now turn our attention to the solutions of the vector field equation and its cosmological 



implications. The equation of motion (10) resembles that of a scalar field, with additional 
contributions to its effective mass stemming from the expansion of the universe. We shall 
study the case in which this mass is dominated by the expansion, ^ V\ which we 
generally expect to happen at early times, and also the case in which the mass is dominated 
by the potential term, <^ V , which we typically expect to happen at late times. We take 
the expansion history as given, which in general requires the vector field to be subdominant. 
In some cases, if the equation of state of the vector field is compatible with the assumed 
expansion of the universe, our solutions also apply to periods in which the vector field is the 
dominant component in the universe. 

In the limit ^ V, the equation of motion of the field takes the form 

fi + 3g. + 3(^-^-;-^^'g-ga±iglg^y,^0. (14) 

The behavior of the solution crucially depends on the size of the term in parenthesis, which 
is typically of order H^. To find explicit solutions it is convenient and realistic to assume 
that the universe expands as a power law, a oc t^, which leads to the exact solution 

M(t) = Cit(^-3p+^)/2^C2t('-'P-^)/', where (15) 
/ = ^(1 + 3p)2 + 1^ [/?i3(p - 1) + 2/34(2p - 1)], (16) 

and Ci and C2 are two integration constants. To be more specific, let us consider for example 
radiation domination, p = 1/2. Keeping the dominant mode only we get 



u 



(t) ^ where q=^2h- 12/3i3//3. (17) 



Thus, the field grows for g > 1, decays for g < 1, and oscillates logarithmically with 
decaying amplitude if q is purely imaginary. Its equation of state can then be calculated by 



substituting (15) into equations (13) and (12), though we shall not do so, as the resulting 



expression is not particularly illuminating. A quite interesting dark energy model in which a 



massless vector field evolves at early times following equation (17) has been recently proposed 
in [101. 
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B. Solutions at intermediate times 



Generally we expect the condition H"^ ^ V to break down at late times, either because 
H steadily decreases, or because the field evolves to regions of larger potential slope as the 



universe expands. In the limit <^ V , the equation of motion (10) resembles that of a 
minimally coupled scalar field, 

V 

u + mil + — u = 0. (18) 

P 

Thus, even if the field is initially frozen, during a small time At it changes by an amount 

It follows that during a Hubble time At = the relative change of the field is of order 
V / {jSH"^). If the latter is large, the field rolls down until the potential slope is sufficiently 
small. 



C. Solutions at late times 



If the potential V has the form ([T]), the equations of motion in the absence of matter 
have fiat space and = (M, 0,0,0) as a solution. We would speak of this configuration as 
the (classical) "vacuum" of the theory. In an expanding universe however, A'^ = (M, 0,0,0) 



does not solve equation (10), unless spacetime is de Sitter and P13 = — 4/?4. To study the 



properties of cosmological solutions around the vacuum, let us substitute into equation (10) 
the ansatz u = M + 5u, with 5u <^ M and V given by equation ([TJ. Keeping only the 
leading terms we get 

/3[Su + 3H6u] - 4XM^6u = 3M[(/?i3 + Af34)H^ - (/? - As - 2(3^)H]. (20) 

The configuration u = M is stable if — 4AM^//9 > 0, which is what we assume in the 
following. In this case, the equation describes a damped harmonic oscillator with a time- 
dependent external force. To solve the equation we shall consider the limit of large masses 
or times, that is H/M ^ 1 (we assume that H ^ H"^.) In this limit, one can think of the 
driving force as being constant, so an approximate particular solution of the inhomogeneous 
equation is 

^«ad = ^ M. (21) 
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This is the leading correction to the solution m = M in powers of the adiabaticity parameter 



H/M. Substituting equation (21) back into (20) one can indeed verify that corrections 
to this solution are suppressed by two additional powers of H/M. The derivatives of this 
adiabatic solution are also suppressed by powers of H/M . Not coinci dentally, in the limit 
H/M <^ 1, the inequality — 4AM^//5 > follows from one of the conditions for classical 



stability at low momenta that we derive in Section VC 

Having found a particular solution of the inhomogeneous equation, we can come back to 
the general solution of the homogeneous equation. 



where B is an arbitrary amplitude and if an arbitrary phase. During a Hubble time the field 



performs y —\/ I3M / H ^ 1 oscillations, while its amplitude remains essentially constant. 



We have seen that the full solution of equation (20) has the form 



U = UsA + SUosc, (23) 



where Suosc is given by equation (22), and, to lowest order. Mad ~ M. In order to find 



the energy density and pressure associated with the full solution (23), we substitute it into 



equations (13) and (12). Because we are interested on cosmological timescales, we consider 
time averages of p and p over Hubble times. Among the terms quadratic in -Uad, those 
without derivatives dominate. Terms that involve mixed terms Mad ^^osc average to zero, 
while the averages of terms quadratic in 6uosc give non-zero contributions proportional to 
the square of the amplitude B. We can thus split the contributions to the energy density in 
an adiabatic component, and an oscillatory one. 

To leading order, the energy density and pressure of the adiabatic solution are. 

Pad = -3(3/5 - 2/5x3 - 2(3,)H'M\ (24) 
Pad = (3/3 - 2/5i3 - 2f3^){3H^ + 2H)M\ (25) 

The subleading terms are suppressed at least by a factor of order H'^/M'^, and, hence, their 
contributions to p and p, of order H^, are negligible. Because the Friedmann equations in 
the presence of matter can be cast as 

3H^ = STrGcPm and - {3H^ + 2H) = SnG^Pm, (26) 

12 



an observer in such a universe trying to determine Gc through the expansion of the universe 
would find a renormahzed value of Newton's constant, 

G 

^ 1 + 87rGM2 (3/5 -2/5i3 -2/54)' ^^"^^ 



As we discuss in Section VC stability requires that 3/5 — 2/5i3 > 0. 



Exactly the same renormalization occurs in aether models We show in Section IV 



that for generic parameter values Gc also differs from the gravitational constant measured 
in "local" solar system measurements Gn- Following the analysis of [37], one can then 
derive a not too stringent order of magnitude limit on the scale M based on nucleosynthesis 
constraints, 

SvrGM^ < 1. (28) 

Since our analysis assumes that M < A, and presumably A is smaller than the scale of 
quantum gravity Mq = (87rG)~^/^, this is also a necessary condition for the validity of our 
approximations. Hence, we typically expect Gc ~ G. Nevertheless, for some parameter 
choices G and Gc can differ significantly. If 8ttGM'^{3(3 — 2/5i3 ~ 2/^4) is fine-tuned to be 
close to —1, the Planck scale one would derive from Gc would be much lower than the actual 
quantum gravity scale Mq- On the other hand, if SttGM^ ^ 1, the actual scale of quantum 
gravity would be much lower than the one guessed from the value of Gc- Of course, this 
would require A ^ Mq- 

The density and pressure of the oscillating solutions are dominated by those terms with 



the maximum number of derivatives acting on u. Thus, substituting equation (22) into 



(13) and (12) and averaging over time we find that the oscillating vector field behaves like 
non-relativistic matter, 

Posc~4AM2— and Posc ~ 0. (29) 

The effective pressure is consistent with the scale factor dependence of the energy density. 
Since we assume that the vector field only couples to gravity, there is no additional decay 
channel through which the vector field could dissipate its energy (if /?i = the spectrum 
of the model consists of a single scalar.) Thus, an oscillating vector field is a potential 
dark matter candidate, somewhat similar to the axion, though we would have to study how 
structure develops to confirm that the vector actually mimics dark matter. These oscillations 
are absent in aether models, because the constraint on the vector field norm only allows the 
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lowest order adiabatic component. Note that in the hmit in which A — oo, for fixed M , the 
adiabatic approximation u^d ~ M becomes exact, and the energy density of the oscillating 
solution diverges, unless 5 = 0. 

The energy density of the oscillating field depends on the value of the amplitude and 
therefore, constraints on the scale M based on this solution depend on initial conditions. 



Nevertheless, based on the adiabatic solution (21), which does not depend on initial condi- 
tions, it appears reasonable to assume that when H ^ M, the amplitude of the oscillations 
was (B/aM) ~ M/X. Under this assumption, demanding that the energy density of the 
oscillating field during nucleosynthesis be no more than the total energy density we arrive 
at the bound 

M < (A'Ml T^^^i)'^' ~ A'/^ 10^° GeV, (30) 
where Tnuci ~ MeV is the temperature during nucleosynthesis, and Mp ^ 10^^ GeV is the 



reduced Planck mass. For A of order one this limit improves the bound (28) by eight orders 
of magnitude. 



IV. NEWTONIAN LIMIT 



The presence of the vector field not only affects gravity on cosmological scales, but also 
locally. In a massless vector field theory for instance, the vector field determines the grav- 
itational constant that enters the Newtonian limit of the theory, and the post-Newtonian 
parameters of the theory generically differ from those of general relativity [36]. In a massive 
theory, one would expect the effects of the vector field to be negligible at sufficiently large 
scales, but it turns out that, even on large scales, the vector field leaves its imprint on the 
effective gravitational constant that appears in the Poisson equation. In this section we 
study some aspects of the Newtonian and Post-Newtonian limits of the unleashed aether. 
A general discussion of the Post-Newtonian limit in the gravitational sector of Kostelecky's 
standard model extension [321 [33] can be found in [28] . 

We would like to determine the gravitational field created by a non-relativistic and static 
source sitting in otherwise empty space. For a potential like ([T| the equations of motion in 
the absence of matter have Minkowski spacetime as a solution, 

ds"^ = -df + dx'^, = (M, 0,0,0). (31) 
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We expect the presence of a sufficiently small and static matter source, with T^,^ = PmS^o^uo, 
to modify this solution slightly, so we can use first order perturbation theory to calculate 
its gravitational field. Because the energy density pm is a scalar under spatial rotations, in 
this limit it suffices to restrict our attention to the scalar sector. Let us hence perturb 

ds'^ = -{l + 2^)dt'^ + {l-2^)dx'^, A" = + {u,dis), (32) 

and solve the equations of motion to ffist order in the perturbations. If we were interested, 
say, in gravitational radiation we would have to study the vector and tensor sectors as well. 
In Fourier space the /x = and n = i linearized field equations ([s]) are 

(/3iA;2 -4AM2)(u + M$) = (/3i3 - /3i - 2/^4) /c^M$ - 4/^4 A;^M(^ - $), (33) 
Pk'^kiS = 0. (34) 



Because we are interested in solutions that vanish at infinity, it follows from equation (34) 



that s = 0. As we discuss in Section VC, stability in this background requires either /3i = 0, 



or /?i > and A < 0. Hence, the relative signs of the operator multiplying u + M$ in 



equation (33) are correct. The same equation can be combined with the linearized i — j 



Einstein equations 

($ - ^)(1 - IGttGM^/S^) = 327rGf3iM{u + M$), (35) 
to solve for the value of Eddington parameter 7, implicitly defined by the relation \E' = 7$. 



If = 0, 7 = 1 follows directly from (35), but in the general case the solution is not that 
simple. 



_ - 4AM^)(1 - IQttGM^P^) - 327TGM^e(34{f3i3 - A + 2/^4) 
^~ (/?iP - 4AM2)(1 - 167rGM2/?4) - USTrGM^k^pj ' ^ ^ 



Inspection of equation (36) reveals that the Eddington parameter is scale-dependent, and 
that the scale in the problem is set by the squared mass of the scalar, — AM^. Therefore, for 
coefficients /5j of order one there are two relevant regimes: k^ <^ |AM^| (large scales), and 
k"^ ^ |AM^| (small scales). The value of 7 in both regimes is 

1 large scales, 

7 = <! P^ + 167^GM^P,i(3,-2(3^^-Ai3,) (37) 
A-167rGM2/34(A + 8/34) ' 
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On short scales the value of 7 agrees with the one calculated for massless vectors in [36] . 
while on large scales it approaches the value in aether models [101 IH] and general relativity 
7 = 1 (the case Pi = = requires special consideration, and is analyzed below.) At 
present, the best measurement of 7 comes from Doppler tracking of the Cassini spacecraft 
[12], which gives 7 — 1 = (2.1 ± 2.3) x 10"^. In the limit A — > cxd, all lengths fall in the long 
wavelength regime. 

It is also instructive to study how u is related to the potential $. Combining again 



equations (33) and (35) we can solve for u in terms of $. To keep the expression manageable, 
let us write down the solution in terms of the linear perturbation of a quantity proportional 
to the perturbation of the squared norm A^A'^, 

A;2[/5i-167rGM2/54(A + 8/54)]-4AM2(l-167rG'M2/54) ' ^ ' 

In the massive case the vector field perturbation u is still non-negligible even far away from 
the matter source. In fact, at long distances, u approaches — M<1>, which typically decays 
as 1/r, where r is the distance to the source. This is quite different from a massive scalar- 
tensor theory, in which to first order the scalar field perturbation 6ip satisfies the equation 
[k"^ +m?)5(p oc pmi and thus vanishes exponentially far away from matter, 5(j) oc exp(— mr)/r. 
Clearly, the origin of the difference is that u + M $ plays the role of the short-ranged scalar 
5 if. 

Once we know how \1/ and u are related to $ we can proceed to solve for the gravitational 
potential $. The linearized 00 Einstein equation is 

= AttG { [(2/?i3 + 4/?4)A;^ - SAM^] M{u + M*) - 2pi3k^M^^ + Ap^k'^M^^ - $) + Pm} 

(39) 



Inserting \1/ = 7$ and equation (38) into (39) leads to the Poisson equation A$ = ATrGcsPm- 
However, the effective gravitational constant Ges is actually not constant, but scale depen- 
dent, 

( [1 - 8ttGM^{Pi + 2/^4)]"^ large scales, 

^ I Pil - 87rGM2[2A(/3i3 - P,{l - 2)) - (A3 + 2/?4)(/?i3 + 2P,{1 - 27))] '""^^^ 

(40) 

The vector field again renormalizes G on large and small scales, but by a different amount 



than in the cosmological case (27), even on large scales. This disagreement in the local 
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and cosmological values of G points to a violation of Birklioff 's theorem in these theories, 



and also is what led to the nucleosynthesis constraint discussed in Section III Note that 
the change in the gravitational constant caused by the term proportional to (3^ is the one 
that follows from replacing — in the action ([2]), and is the same both for the 

cosmological and Newtonian gravitational constants. If we set (S^ = 0, our value of Gefr on 
large scales agrees with the one found for aether models in j37]. Note that the value of f3 
does not have any impact on the Newtonian limit of the theory. 

As an illustration of these results, let us consider a massive vector with = 0. It follows 



automatically from equation (37) that 7 = 1, so in this case the theory agrees with the 
current measurements of 7 both on short and large scales. The effective Newton constant 
takes the simpler form 

Gcs I [1 - 87rGM2/5i]-i large scales, ^^^^ 

G I [1 - 8nGM'^Pis{2 - pis/Pi)]-^ small scales. 

This equation is remarkable because it says that the strength of gravity changes with scale. 
It is hence tempting to speculate a connection with dark matter or even dark energy. For 
instance, if /3i//5i3 > gravity is stronger at large scales than at small scales, so if the vector 
had a Compton wavelength of a few kiloparsecs this could eventually explain the flatness of 
galactic rotation curves without dark matter. Recall that if /?i 7^ 0, stability imposes (3i > 
and /5i3 > 0. 

Such a scale-dependent behavior of the gravitational coupling does not occur in massless 
vector field theories or aether models, though it does happen in theories with more com- 
plicated Lagrangians |13]. In connection with modified Newtonian dynamics a specific 
model that involves a fixed-norm vector has been proposed by Bekenstein in Ideas 
about degravitation on large scales have been explored in |36] and references therein. 



Although our expressions (37) and (40) still hold when jSi = or jS^ = 0, they become 



ill-defined when both /3i and (3^ vanish. This case however is particularly simple. It follows 



directly from equation (35) that 7 = 1, while the remaining linearized equations lead to 

GeS 



4AM2 



G 4AM2 - fc2 



(42) 
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Therefore, on short scales the theory does not have a Newtonian hmit (the gravitational 
potential satisfies the equation oc p), while if the field is massive, G^s ~ G on long 
scales. 



V. GHOSTS AND INSTABILITIES 

The strongest constraints on the parameters in the action ^ are set by the requirement 
that the theory be classically and quantum-mechanically stable. In this somewhat technical 
section we proceed to analyze the conditions that classical and quantum-mechanical sta- 
bility place on the parameters of the action (|2]). Analogous conditions for vector fields of 
constrained norm have been derived in |17l HE] , while a somewhat similar analysis of mass- 
less vector fields of arbitrary norm has been carried out in [49j. Instabilities in models of 
infiation supported by vector fields have been recently reported in [501 EB IS2] • Our stability 
analysis closely resembles the one for massive gravity in |53] 

In order to analyze the stability of unleashed aether models, we consider vector field 
fiuctuations &^ around the background value of the field, 

Af, = A^ + b^. (43) 

Because we are mostly interested in cosmological solutions we assume that the background 
vector A^^ has the form ([o]), and leave the case of a field that points along a spatial direction to 
Appendix |Xj We neglect metric perturbations, which relies on the assumption that A <^ Mq- 
In conformal time coordinates, the metric and background field are hence given by 

ds^ = a^{r]){-drf + dx^), = a"^ ■ (A, 0, 0, 0). (44) 



Note that A plays the role of u in Section III, and is hence dynamical. 



To obtain the Lagrangian of the fiuctuations we insert the expansion (43) and the 
background (44) into the action ([2]). Expanding to quadratic order in the fiuctuations we 
arrive at 

A = - l3id^Kd%'' -{(3- I3{){d^h^f - 4(3{n%l - nbod^b^) + 

+ PisRt.ub^b'' + a^p^Rh^h^" - a^Vb^b^" - 2a^V" ■ (v4^&^)^ (45) 

where indices are raised with the Minkowski metric rj^^ , Ti = dlog a/di], and a prime denotes 
a derivative with respect to A^A'^. For arbitrary parameters (3i and (3, the four components 
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of the vector field are dynamical, because the determinant of the Hessian matrix 



(46) 



d[dh^/d'r])d{dhy/dri) 

is non-zero [54J. However, as we shall see below, for /5 = or /5i = the determinant of 
the Hessian vanishes, signaling that some of the vector field components are constrained. 
This is what happens for instance for the gauge invariant Maxwell theory (/5 = 0). We thus 
study the degenerate cases /3 = and /5i = separately. In aether models, the fixed norm 
constrain requires 6° = [17j, and thus eliminates one of the four dynamical fields. Note by 
the way that if we expand the potential of equation ([T]) to quartic order in the perturbations 
we obtain an additional term 

D Xih.h'^f- (47) 
Thus, for A > 1, the theory is strongly coupled. 



A. Stability in the vector sector 

The background value of the vector field determines the set of unbroken spacetime sym- 



metries. Because points along the time direction, the action (45) is symmetric under 
spatial translations and rotations. It is convenient to decompose the fields in irreducible 
representations of this unbroken symmetry group. Expanding as usual 

ho = —u and hi = diV + Vi, with diV^ = 0, (48) 

decouples the two scalar modes u and r from the two vector modes in v. The Lagrangian of 
the vector sector is then given by 

£^ = -(3^df,v ■ d'^v - a^ml v ■ v, (49) 

where 

ml = V'- 3(/3i3 + 4(3^)H' - (A3 + QP4)H. (50) 



The factor of a in the mass term of equation (49), as well as the contributions to the mass 
proportional to the curvature explicitly break time translation invariance. Therefore, the 
quantity one would naively identify as the frequency of the modes in the vector sector 

=P^ + -w^nl (51) 
Pi 
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Vector 


Low p 


High p 


Classical 


Piml > 


/ 


Quantum 


/3i >0 


Pi>0 



TABLE L Stability conditions in the vector sector. The checkmark means that the condition is 



automatically satisfied. The mass rriv is defined in equation (50). 



also depends on time (in this regard one can consider a^m^//?i as the "mass" of these 
excitations.) We can nevertheless quantize the theory under the approximation that the 



frequency is constant if the relative change in u during the characteristic time Ar] 



IS 



small 



1 dio 
oj'^ drj 



< 1. 



(52) 



In the short wavelength regime, •p'^ /a? ^ m^, this adiabaticity condition is automatically 
satisfied, but in the long wavelength regime p'^/a^ <^ m^, it imposes restrictions on the 
magnitude of V and its time derivatives. For instance, if V <^ H^, the mass is of the 
order of H^, and hence uj~'^{dujv/ drj) is of order one. Our analysis in this section exclusively 



applies in the adiabatic regime (52). In this regime, we can set all explicit factors of a to 
one. 



Returning now to equation (49), and expressing the vector in terms of two orthogonal 
transverse polarizations immediately reveals that the vector sector is ghost free for /5i > 0, 
and free of classical instabilities at low momenta (p^ <^ I'^^l) for ^^/A ^ 0. These condi- 
tions are summarized in Table |T} We should also point out that tachyonic, long-wavelength 
instabilities in a theory may be allowed, as long as they are not too severe. Their severity 



depends on the magnitude of the tachyonic "mass" uj"^ 



-m 



which determines the rate 



of growth of the mode, proportional to exp(mj?7). Since we mostly assume that all masses 
are much larger than the Hubble scale, the presence of a tachyon would imply instabilities 
on timescales much shorter than the Hubble time, which is likely to be problematic. On 
the other hand, tachyonic masses of the order of H (or smaller) should not be much of a 
problem. Note that if /5i = 0, the vector sector is not dynamical. 
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B. Stability in the scalar sector 



To proceed further, we expand the scalar fields in appropriately normalized Fourier modes, 



u 



u{p) exp{ipf^x^) 



d'^p r{p) 



exp{ip^x^ 



(53) 



(27r)4 -v-^---v-A^- ^' J (27r)4 \p\ 

Since the fields are real, it follows that u*{p) = u{—p) and r*{p) = r{—p). Substituting the 



expansions (53) in (45) and assuming as above that all explicitly time-dependent quantities 



are constant we arrive at the action (in Fourier space) 



5. 



d^p 




(54) 



where the matrix D is given by 

I /3u^-f3ip^ + ml 
\zuj\p\{f3-(3,)-2(3\p\H 

and we have also defined 



-iu\p\{P-P^)-2p\p\H 
-I3i uj^ + (3p^ + ml 



(55) 



ml = 2V"A^ -V' + {213 + 3/3i3 + l2p^)H^ - (2/3 - 3/3i3 - QPa)H. 



(56) 



In spite of our notation, the mass parameter is not the mass of the field u, because 
the Lagrangian is not diagonal in that field. Note that in an expanding universe, the terms 
proportional to R^^ and R in equation ^ also contribute to m„ and m^. These contributions 
are typically negligible, unless the derivatives of the potential are sufficiently small. 

The inverse of the matrix D is the field propagator. Thus, in order to find the propa- 
gating modes we just have to find the values of cu^ at which its eigenvalues have poles, or, 
equivalently, the values of uj"^ at which the eigenvalues of D have zeros. The determinant of 
D is 

detD = -(5f3iij^+{2pf5if-f3iml+Pml)u^-{f3f3if-(3mlp^+pimlp^+Ap^^ 

(57) 

Requiring that det D vanish we thus arrive at the frequencies of the two propagating scalar 
modes 



2 ^2 I 

UJ = p + 



"""iA, where 



2/?i 2/? 



A 



^ 4 V /? ^ A 



(58) 
(59) 



21 



Observe that some of the factors of iJ^ in the determinant (57) are multiphed by the mass 



terms and m^. Since factors of lo arise from time derivatives, and because the mass 
terms contain derivatives of the potential, we may say that the potential is also part of the 
vector field kinetic terms. 

The theory is free of classical instabilities (exponentially growing modes) if the frequency 



uj in equation (58) is real. In the limit of low momenta, in which p is negligible, the two 



frequencies approach^ 



Co- , 



0J_ 



2 

/3i 



+ 



+ 



1 - 

1 + 



(/?-/ji)(m^ + mg)-4/j^g^ 

/3im2 + /3m2 
(/5-A)(m2+m2)- 4/3^/721 



p 2 + O {p^) 
p' + O {p') 



(60) 
(61) 



in which, for later convenience, we have included terms that are typically subdominant at 
low momenta. In the limit of high momenta, in which mass terms can be neglected, the 
frequencies become 



p^ ± \p\ 




1 



{ml + ml) 



4- 



(62) 



Demanding that the mode frequencies be real at low and high momenta we arrive at the 
conditions on Table[Tl| which if satisfied actually imply that the system is classically stable for 



all momenta (this follows from a detailed analysis of equations (58) and (59).) Inspection 



of equations (58) and (59) shows that the regimes of high and low momentum not only 



depend on the masses of the fields, but also on the values of (3 and jSi. Note that in the 
high-momentum limit, the leading term in the frequency gives uj± p"^. We consider the 
subleading contributions, proportional to because if the latter are imaginary there are 
instabilities on arbitrarily short timescales. 

The residues at the poles of the eigenvalues of determine whether the propagating 
modes are ghosts. If these residues have negative signs, the corresponding modes have 
negative energies. In the presence of ghosts, the vacuum can then decay into positive energy 
quanta and negative energy ghosts, leading to a quantum-mechanically unstable theory. How 
severe this instability is depends on the way the ghosts couple to positive energy quanta, 
and on the volume of phase space available for the decay [SHI EZ] • The coupling of ghosts to 



We are assuming that (m^//3 + ml//3i) < 0. Otherwise the vahies of uj+ and ti>_ become interchanged. 
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positive energy quanta is essentially determined by the coupling constants of the theory and 
the magnitude of the residues, while the phase space available for decay typically depends on 
the cut-off up to which the theory is valid. Choosing any of these parameters appropriately 
might tame ghost instabilities, although we shall not explore this possibility here. In this 
context, it is worthwhile pointing out however that with Lorentz-invariance spontaneously 
broken, the cutoff of the theory involves spatial momenta rather than Lorentz-covariant four- 
momenta. To illustrate this let us consider one of the many higher-dimensional operators 
that we neglected in ([2]), 

^AMM''A-(V^V,AO(V,V.A,), (63) 



where (5 is a dimensionless parameter of order one. Expanding (63) to quadratic order in 
the fluctuations 6^ in Minkowski space we obtain 

D S^u\ (64) 

Hence, this term would make an effective contribution to the matrix D which would become 
comparable to the lowest order term (3uj'^v?' if 

« ^. (65) 

As we have seen, at low momenta the dispersion has the form uj"^ = m^, where, by definition, 
m is the mass of the field, while at high momenta the dispersion relation typically approaches 
uj"^ = c^sP^i where is a coefficient of order one. Thus, our effective description is certainly 
not valid for masses and spatial momenta beyond / A^. In fact, higher-dimensional oper- 



ators analogous to (63), of the form 

_L_A/^i . . . A^-A^^ ■ ■ ■ A^"iV,, ■ ■ ■ V,^A-){V,, ■ ■ ■ V,^A^), (66) 

point to a breakdown scale of the order A'^/A. Note that in the limit A ^ 0, the corrections 
due to the new term goes to zero. In that limit Lorentz invariance is restored, and large 
spatial momenta can be made small by a Lorentz transformation.^ 



^ A similar phenomenon was described in the context of a scalar field theory in an inflating universe in |58j . 
There, the scale of Lorentz symmetry breaking is set by the Hubble constant, which plays the role of A 
here. 
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If A+(ct;^) and A_(ct;^) are the two eigenvalues of the matrix D, and the former have zeros 
at, respectively, a;+ and the residues at the poles of D^"^ are^ 



lim 



1 



(67) 



The theory is free of ghosts if both residues are positive. Since the two eigenvalues of the 
matrix D are 



Al,: 



iiD ± Vtr^D -4detD 



the condition for the absence of ghosts is then 



1 



1 ddetD 
t^D duj"^ 



> 0. 



(68) 



(69) 



The trace of the matrix D at uj± is 
trD\.=2iP-P,)p' + 



whereas the derivative of the determinant is 

ddetD 



^+-)±(/^-A)A, 



T2/3/?iA. 



(70) 



(71) 



From the previous equations it immediately follows that the theory always contains a high- 
momentum ghost. Indeed, in the limit of high momenta, the trace approximately equals 
the same value both at uj+ and uj-., while the value of ddetD/duj'^ changes sign. Therefore, 
the residues cannot be positive for both frequencies. More precisely, in the limit of high- 
momentum the residues approach 



(72) 



which implies that one linear combination of fields is a ghost, while the remaining orthogonal 
combination is stable (classical stability forces the term inside the square root to be positive.) 
Because at low momenta the residues approach 



1 



Pi, 



(73) 



® In some cases, one of the eigenvalues may have multiple zeros. In that case, evaluation of equation (67 1 
at those zeroes gives the corresponding residues. 
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Scalar / 0, /3 / 0) 


Low p 


High p 


Classical 


Pi-ml>0, i3-ml<0 


(/Jfi - p-')iml + ml) - 4(3I3^'H^ > 


Quantum 


^1 > 0, /3 < 


# 



TABLE IL Stability conditions in the scalar sector for (3i ^ and /3 7^ 0. The cross means that the 
condition cannot be satisfied. The parameters mv and are defined, respectively, in equations 



(50) and (56|. 



it is possible to avoid the existence of low-momentum ghosts for appropriate choices of f3 
and Pi. Classical and quantum stability conditions are summarized in Table [IT] 

Note by the way that there is a ghost at high momenta for f] = (3i, even though equation 



(72) breaks down for those values. In this case, the trace of D has always the same value, 



irrespectively of p. If H or any of the two scalar masses is non-zero, A is non-zero, and, as 



before, equation (69) cannot be satisfied. This is immediately apparent in the Lagrangian 



(45). For (3 = Pi, the Lagrangian is already diagonal in the fields b^, and because of the 



indefinite signature of the Minkowski metric, at least one of the fields has to be a ghost. 



1. Massless case 



It is important to realize that the high-momentum limit is different from the strict mass- 
less case 



0. 



(74) 



Indeed, in the massive case there are always two propagating modes in the scalar sec- 
tor at high momenta, with frequencies that approach the relativistic dispersion relations 



p^{l ± 0{m/\'p\)). In contrast, in the massless case, A in equation (59) vanishes, and 



equation ( 71 ) implies that there is a double pole at oj"^ = p ^ instead of two single poles at 



UJ 



p"^ . Thus, in the massless case there is a scalar with higher time derivatives plus the 



two degrees of freedom of the vector sector. Lagrangians with higher time derivatives are 
usually quantized using a procedure developed by Ostrogradski, and are typically unstable 
(see [5^ for a review.) 

The same conclusion follows by substituting the Lorentz-covariant decomposition 
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= d^(f) + Vfj,, with dfj.v'^ = into the action 

Sm = j -^F^^F^- -l3{V^A^^f-m^A^A^^ . (75) 

This yields a Lagrangian which describes massive electrodynamics in "Lorentz gauge" and 
a scalar cf) whose kinetic term contains higher order derivatives, 

= -^{d^v, - d,v^){d^v'' - d'v'^) - /?(n0)2 - m'd^cpdy - m\^v^. (76) 

In this form, it is manifest that the Lagrangian not only describes a vector field, but also 
a scalar field (under Lorentz transformations.) By inspecting the scalar kinetic term in 



equation (76) one would naively conclude that in the massive case there are two scalar 
poles, at J? = —im?/(3 and = 0, but this is illusory. Because df^v^ = 0, the propagator of 
the original vector field is 

PiP + m"' \ P / pp^ + p^ 

The propagator of v is just that of electrodynamics in gauge [6^, with ^ — and 



an added mass term. It can be obtained by inverting the kinetic term in equation (75) 
and letting (3, which plays the role of 1/^, tend to infinity. This propagator is manifestly 
transverse (projects onto the subspace perpendicular to p^), while the contribution from </> 



to A A is clearly longitudinal (projects onto the subspace parallel to p^.) In the form (77) 
it is now clear that there are no poles at p"^ = 0. Instead, there is a pole at p"^ = —im?/j3 
in the longitudinal sector, and three poles at p^ = —im?/(3i in the transverse sector. These 
correspond to the same frequencies that we found in our previous analysis, as the reader 



can readily verify by setting H = 0, = m and = — m in equations (58) and (59) 



Note that we could have also obtained ( 77 ) by directly inverting the differential operator in 



equation (76). 

Let us finally discuss = 0, as we intended to do. In this case, the poles at = that 
canceled in the massive case remain, and give rise to factors proportional to There 
is then a double pole at = 0, and two single poles at = 0. These considerations also 
illustrate that a free (V^ = 0) vector field in fiat space can behave very differently from a 
free vector field in an expanding universe if /3, /3i3 or (3^ are different from zero. 
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2. Degenerate Limits 



Our conclusions do not apply to the two degenerate limits in which equation (58) breaks 
down namely, for /5 = and (5\ = 0. In both cases, the determinant of the Hessian matrix 



(46) vanishes, or equivalently, the determinant of D, equation (57), is not quartic in u, but 
quadratic. The Lagrangian therefore describes only one scalar propagating mode, instead 
of two. 

When (3 = 0, the "kinetic" part of the Lagrangian is that of the electromagnetism and 
the Proca Lagrangian. The determinant of the matrix D only has one zero at 



UJn 



ml 



mt 



(78) 



and the residue at this pole is 



+ 



2 2 

m.. — mt 



m: 



■P' 



(79) 



Thus, equations (78) and (79) result in the stability conditions in Table III Note that some of 
these conditions are the same as in the vector sector. Table |T| and that if the theory is stable 
both in the low and high momentum limits, it will be stable for all momenta. If = 0, 
the equation det D = has no solution and the theory does not have any propagating scalar 
degrees of freedom. From the point of view of flat space, equation (77), setting /3 = 
eliminates the longitudinal polarization associated with the scalar (p. In the massless case, 
one of the the remaining transverse polarizations disappears, leaving just the two modes of 
the vector sector. In this case, the theory has an additional symmetry, gauge invariance. 

As illustration, let us consider the familiar example of the Proca Lagrangian, in which 
/?! = 1/2, /? = /5i3 = /?4 = and V = m^A^A^/2. For a vanishing background field this 



implies ml = m /2 and m^ = —m /2. It follows from Table III, that the Proca field is stable 
for > and both classically and quantum-mechanically unstable for m^ < 0. Again, 
the massless limit is not continuous: for m^ = there are no propagating scalars. Massive 
versions of electromagnetism with mass terms that violate Lorentz invariance (m^ 7^ ~^V) 
have been studied in [6II |62lj . 

If /3i = the determinant of D vanishes at 



mz 



m: 



p , 



mt 



(80) 
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Scalar {(3 = 0) 


Low p 


High p 


Classical 


Piml > 


mlml < 


Quantum 


/3i >0 


ml-ml> 



Scalar = 0) 


Low p 


High p 


Classical 


f3ml<0 


ml{ml-4l3H'^) < 


Quantum 


p <0 


ml -ml- 4pH^ < 



TABLE HL Stability conditions in the scalar sector for timelike background fields in the degenerate 
cases /3i = or /3 = 0. If the theory is stable both in the low and high momentum limits, it follows 



from equations (78), (79), (80) and (81 ) that it is actually stable for all momenta. The parameters 



m^ and are defined, respectively, in equations (50) and (56). 
while the residues at the poles are 



7 _ 1 ml 



ml 



(81) 



Combining both (80) and (81) we arrive at the stability conditions in Table III If these 
conditions are satisfied, the theory is stable for all momenta. Once again, if ml = the 
inverse of D has no poles, and the theory does not contain any dynamical scalar (and, since 
/3i = 0, there is no dynamical vector either.) In a Lorentz-covariant language, setting (3i = 



eliminates the three transverse poles in the propagator of equation (77). In the massless 



case, the longitudinal polarization also disappears, leaving a theory without dynamics. In 
fact, in this limit the theory has an additional symmetry —>■ + ^^I'^'^^d^Cp^, where Cp„ is 
an antisymmetric tensor (a two-form), which can be used to gauge all the components of 6^ 
away. 



C. Stability around the minimum of the potential 



As an important application of this analysis, let us consider stability in the case where 
the vector field sits at the minimum of a potential given by equation ([T]). As discussed 
in Subsection |III C this is generically possible only in fiat space. Substituting H = and 
equation ([T| into equations (50) and (56) we find 



and ml = 4XM' 



(82) 



Let us consider the case /3i 7^ and /? 7^ first. The spectrum then contains two vector and 
two scalar degrees of freedom, with dispersion relations at low momenta given, respectively, 
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by equations (51) and (60), 



2 -'2 2 -'2 2 4AiVf 2 ^2 lor>\ 

^v^=P , ^V2=P ^ ~ ^s^ ~ JP ■ (83) 

The three "massless" excitations fi,f2 and s_ can thus be interpreted as the Nambu- 
Goldstone bosons associated with the spontaneous breaking of the three generators of 



Lorentz boosts. At low momenta, the corresponding residues are, from equations (49) and 



(73) 



Z,, = Z,, = Z,^ = -1/13, Zs_ = l/Pi. (84) 

In consequence, at low momenta the theory is either classically or quantum-mechanically 
stable, but not both. The absence of ghosts requires Pi > and P < 0, which then implies 
that s_ is classically unstable. Conversely, classical stability demands A//5 < and P/Pi > 0, 
which then implies that either vi,V2 and s_, or s+, are ghosts. The only way to escape this 
conclusion is to choose a well-behaved set of massless modes (/5 > and Pi > 0) and push 
the massive mode — either a tachyon or a ghost — out of the domain of validity of the effective 
theory. Indeed, as mentioned above, if the mass of is of the order A^/M, we cannot trust 
our analysis anyway, so the presence of a tachyon or ghost in the spectrum is not necessarily 
fatal. In principle, this pathology could be solved by the ultraviolet completion of the 
theory. Note that these instabilities cannot be avoided by lowering the mass scale M and 
"squeezing" the low-momentum regime out of the spectrum of the theory; at high momenta, 
^ 4AAf^, according to Table [ll| and equation (72), either s+ or s_ are ghost-like too. 



Formally, we can decouple the massive mode s+ by sending |A| to infinity, which is the 



limit in which we expect to recover aether models. In fact, the low momentum spectrum (83) 



and residues (84) of the massless modes vi,V2 and s_ agree with those of aether theories 
In that sense, unleashed aether models can be regarded as a high-momentum completions 
of the former, which, as we have argued, are inconsistent for generic values of P and Pi, 
unless A is larger than (A/M)^. This is not just a fine-tuning problem. Because we expect 
A/M ^ 1, it also means that the theory is strongly coupled. The only way in which strong 
coupling issues could be circumvented (without avoiding fine-tuning) in a sensible low energy 
effective theory is by considering < /? < (M/A)'*. In any case, from this perspective it 
is also easy to understand the agreement on large scales of the Newtonian limits of both 
constrained and unleashed aether models. At low momenta, both are essentially the same 
theory. 
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These instabilities can be averted by setting (3i = or f3 = 0. Since no non-trivial 
cosmological solutions with (3 = exist, let us hence consider j3i = next. In this case 
the vector sector is non-dynamical, and because ml = 0, there are no dynamical scalars 
in the theory either. This is however a degenerate limit. In an expanding universe, the 



field is slightly displaced from the minimum, as discussed in Section III Assuming that 



^ 4AM^ and inserting equation (21) into equation (50) we find, to leading order 



m2 = -(2/5i3-3/3)if. (85) 

Hence, ml is of order and generally different from zero. The spectrum of the theory then 
consists of a single scalar sq with dispersion relation and residues given by, respectively. 



^ — 1 - — — ~ and Zo^ -- + — ^. (86) 

Thus, the scalar is stable at low momenta if /? < and A > 0. Because H < whenever 
the null energy condition^ is satisfied, the scalar is stable at high momenta if in addition 



2/3i3 — 3/5 < 0. Comparing with the general case above, equation (83), we see that at low 
momenta the effect of setting /5i = is to remove the (classically unstable) scalar Goldstone 
mode S- from the spectrum. The remaining scalar is massive, and thus is absent from the 
low- momentum description provided by aether theories. In fact, aether models with Pi = 
do not have any propagating degrees of freedom. 



The case /3 = is somewhat similar. Using equation (85) we find that the dispersion 
relation and residues of the only scalar mode are 

..^ « -2ft,// (1 - ^) . and Zo«-i^-^ (87) 

Thus, using H < again, stability at low momenta requires (3i > and (3i3 > 0, while 
stability at high momenta then imposes A < 0. 



D. Relation between energy conditions, instabilities and superluminality 

A violation of the standard energy conditions is often associated with superluminal prop- 
agation and instabilities in a theory |56| 163] . In k-field models for example, the absence of 



p > and p + p>0- 
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ghosts implies the null energy condition p + p > 0, and also places at the same time restric- 
tions on the value of the speed of sound (see for instance [M])- It is therefore instructive to 
check whether these relations also extend to unleashed aether models. 



From equation (58) one could immediately calculate, say, the group velocity of the per- 
turbations. 



For instance, it follows from equation (83) that the propagation of s_ is subluminal for 



(3/(3i < 1, and superluminal for > 1. As discussed in Section VC, the stability 



of that mode (in our restricted sense) only requires (3 > and /5i > 0. Similarly, from 



equation (62), at high momenta the dispersion relations of the scalar modes approach the 
relativistic expressions uj± = p, whilst our stability analysis indicates that in that limit there 
is always a ghost in the spectrum. Thus, subluminal propagation seems to be independent 
from stability considerations (in this respect, the unleashed aether is analogous to k-fields.) 
The requirement of subluminal propagation is sometimes used to constrain the parameters 
in this type of theories, but we believe that there is nothing inconsistent with faster than 
light propagation in this background. The constant vector field defines a preferred frame 
and explicitly breaks Lorentz invariance. Even if a signal propagates superluminally, it will 
move forward in our preferred time coordinate, it will not be possible to construct closed 
timelike curves, and no conflict with causality will arise. For different viewpoints on the 
issue see e.g. [651 ESI E?! [68]. Readers concerned about superluminal propagation should 
impose appropriate conditions on the parameters of the theory. For instance, the additional 
constraint < 1 implies that the mode s_ propagates subluminally (at low momenta). 
Analogous conditions in the gravitational wave sector have been discussed in [17j. Note that 
signal propagation in the vector sector is never superluminal. 



In the limit in which the expansion of the universe is negligible the energy density ( 13 ) is 



p = -pix^ + V. (89) 



In the same limit, the sum of pressure and energy density is 

/5l3 + 2/^4 1 \ T// 2 



p + p = 2 l ^'-" - Ij V'u' - 2(/3i3 + mu'. (90) 

Therefore, the energy density is positive for (3 < Q and > 0, while p + p is positive if 
/3i3 + 2/^4 < and V'{l3i^ + 2/^4 - /9)//5 > 0. Examination of Tables |l| to |lTl| shows that 
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these conditions bear a marginal resemblance with the conditions for the absence of ghosts 
in the limit of low momentum, which is the one that applies for an homogeneous field. In 
particular, the absence of ghosts implies /5 < 0, which is the condition for the "kinetic" part 
of the energy density to be positive. Otherwise, at least in this particular case, it appears 
that the weak and null energy conditions have little to do with the presence or absence of 
ghosts, at least in this class of models. 

VI. SUMMARY AND CONCLUSIONS 

At sufficiently long timescales and large distances, any generally covariant theory that 
contains a real vector field coupled to gravity, and is invariant under —* — A^, should 
be described by an action of the form (|2]). Massless vector-tensor theories and aether 
models [12] can be regarded as particular limits of this "unleashed aether" class. We have 
identified some of their common features, and pointed out in what other aspects they differ 
significantly. 

Stability imposes severe conditions on the parameters of the action ([2]). The spectrum 
of the theory generally consists of a massive excitation and three massless fields, which can 
be interpreted as the Nambu-Goldstone bosons associated with the spontaneous symmetry 
breaking of invariance under boosts. At least one of these four fields is always a tachyon 
or a ghost. For appropriate choices of parameters however, it is possible to decouple the 
dangerous modes. Setting (3i = eliminates all the transverse components (in a Lorentz 
covariant way) of the vector. And setting j3 = eliminates just the longitudinal polarization, 
which is dynamical only for a massive vector or in an expanding universe. Both choices result 
in a well-behaved theory. Non-trivial cosmo logical solutions exist only if /5i = 0, that is, if 
the kinetic term of the vector field is proportional to the squared divergence of the vector. 
We have also extended our stability analysis to models in which the vector field points along 
a spatial direction (see the Appendix |A]) and found that they also suffer generically from 
ghosts. 

There is yet another alternative to eliminate the unstable fields from the theory. For 
appropriate choices of parameters, the only unstable field in the theory is massive. If its 
mass is sufficiently heavy, its excitations are beyond the regime of validity of our low-energy 
description and can be ignored. The remaining three fields, the Nambu-Goldstone bosons. 
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are the only fields present in the spectrum of aether models. In that sense the aether can 
be regarded as a low-momentum description of unleashed aether models. 

Since the symmetries of the universe only encompass spatial translations and rotations, it 
is possible for a vector field to have a non-vanishing expectation value along its time direction. 
In unleashed aether models such a non-vanishing cosmic vector field leads to a host of 
interesting gravitational phenomena. In the Newtonian limit of the theory Newton's constant 
is typically scale-dependent, and for some parameter values the Eddington parameter 7 
differs from its value in general relativity. On cosmological scales, the vector field leaves an 
imprint on the expansion of the universe, even in the vicinity of a potential minimum. The 
field renormalizcs Newton's constant, like in aether models, but it also typically oscillates 
around it. Under relatively mild assumptions, these oscillations strongly constrain the value 
of the field at the minimum of the potential. 

Perhaps the most important conclusion is that self-consistent and phenomenologically 
acceptable modifications of gravity are difficult to come by. In this article we have added a 
massive vector field of unknown origin to the gravitational sector, and we have employed a 
low-energy effective approach to study the implications of this addition. We have seen that 
general vector field theories have a very rich gravitational phenomenology, but that self- 
consistency and agreement with basic cosmological and solar system observations impose 
severe constraints on this class of models. Nevertheless, a large volume of parameter space 
is still allowed by our requirements, and many interesting phenomenological aspects remain 
to be investigated. Unleashed aether models can wander far away. 
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APPENDIX A: STABILITY OF SPACELIKE VECTORS 



The stability analysis in Section |V] assumes that the vector field points in the time direc- 
tion, as required by isotropy. But even if the field is aligned along a spatial direction it is 
possible in some cases to restore the isotropy of spacetime by either considering appropriate 
configurations of multiple vector fields [H |12.J, or by simply assuming that the vector is 
subdominant. Our stability analysis can also be extended to this type of configurations, for 
which the background value of A'^ is spacelike. 

If the vector field points along a spatial direction, we can always pick our coordinates so 
that A^^ points in the z direction, 

Af" = a'\0,0,0,A). (Al) 

In order to determine whether such a configuration is stable, we consider again the fiuctu- 
ations around the background value of the field, = + 6^. But instead of expanding 



the fiuctuations as in (48) we decompose the perturbations in irreducible representations of 



5*0(2), the group of rotations that leaves A^ in equation (Al) invariant. 



bo = —u, ha = daT + Va (with So-f " = 0), = t, (A2) 

where a runs over the transverse spatial coordinates, a = x,y. Thus, under S0{2), u, r and 
t are scalars, while v is a vector. 



Substituting the expansion (A2) into equation (45) we obtain the Lagrangian for the 



perturbations. The Lagrangian for the transverse vector modes t;" takes the simple form 

= -Pid^vad^^r - mlvar, (A3) 

where ml is given by equation (50 ) (recall that we assume that a?m1 is adiabatically constant, 
and thus set a = 1.) Thus, the absence of instabilities leads again to the conditions in Table 
[Tj For spacelike background fields the scalar sector in this case is slightly more complicated, 
because there are three coupled scalar fields, instead of two. In Fourier space, their action 
is 

Ss = -\j {[P^'- PiP' + ^l] u*u + [-Piu'p^' + Ap^pI + (/5 - A)p1 + pK\ f*f + 

+ [-/?icu2 + 2 + _ /5,)pJ + rhl] n - ApHp^ru + (A4) 
—Ai(5Hp\\^ i*u — 2{(3 — jdi) [iujp_if*u — upufu — ip±p\\ t*r] + /i.e. |, 
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where we have defined 



^2 



-V + (2/3 + 3/3i3 + 12I3^)H^ - {2(3 - 3/3i3 - Q(3^)H, 



mi = V' + 2V"A' - 3(^3 + ^I3aW ' (As + QI3^)H, 



(A5) 
(A6) 



and p|| = Pz, Rl = (Pa',P?/,0). Note that differs from the in equation (56). For later 
convenience, we have not normahzed s canonically. 

A general analysis of the stability of this system would be rather cumbersome, because 
we would have to diagonalize the 3x3 matrix. For our purposes however it suffices to focus 
on particular values of the momentum p, for which some of the scalars decouple from the 



others. Specifically, setting pj_ = in equation (A4) we arrive at the action 



+ 2i(3Hp\\ [u*t-fu] + {(3 - f3i)ujp\\ [u*t + t*u] }. 



(AT) 



This has the same form as equation ( 54 ) . The field r is non-dynamical and has disappeared 



from the Lagrangian, while t plays the role of r in equation (54). The only difference is 



an overall factor of i in the non-diagonal terms, which does not affect the trace nor the 
determinant of the matrix D. Thus, the stability of perturbations that only depend on z 
(and time) is determined by the same conditions we derived in the previous section, with 



mass parameters given by equations (A5) and (A6). It follows for instance that for generic 



values of the theory has a high-momentum ghost. 

The presence of a ghost at high momenta is not exclusive of z dependent perturbations, 
and also appears for perturbations that only depend on x and y. Indeed, setting py = into 



the action (A4) gives 



S 



(2vr) 



[pu}"^ - (3ipA^ + ml] u*u + [-(3iuj'^p± + I3p^ + p±^m^] f*? 
-2(3Hp±^ [u*r + f*u] - iu{(3 - (3i)p±{r*u - u*r) + 



(A8) 



The field t has decoupled from the other two scalars, whose action has again the form of 



(54) (after canonically normalizing f.) Thus, as before, the theory has ghosts at sufficiently 
high values of p]_ for generic values of /?, Pi. 
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Scalar (/3 = 0) 


Low p 


High p 


Classical 


(iifhj > 0, Piml > 


rh^rhl < 0, m^m^ < 


Quantum 


/5i > 


/?! > 0, m'f - ml > 0, ml - ml > 



TABLE IV: Stability conditions in the scalar sector for spacelike background fields and /? = 0. 



Scalar = 0) 


Low p 


High p 


Classical 


13ml < 


ml{ml - 4PH^) < 0, ml{ml - 4PH^) < 


Quantum 


/3 < 


ml - mj - Aj3H'^ >0,ml-ml- AjdH'^ > 00 



TABLE V: Stability conditions in the scalar sector for spacelike background fields and (3i = 0. 
For (3 = OT j3i = these instabilities disappear if appropriate conditions are satisfied. 



These are listed in Tables IV and |V} As an application, let us consider for instance La- 
grangians with /3 = 0, as in the "new inflation" vector models of or as in the dark energy 



model of [1]. Using equations (50), (A5) and (A6), and neglecting curvature terms we find 



^2 



m; - mz ^ 2V' + 2V"A^ 



ml - ml ^ 2V' 



v. 



(A9) 
(AlO) 
(All) 



It follows automatically from our analysis that all models with "run-away" potentials with 
V < 0, in which the derivatives of the potential dominate over curvature terms, contain 
ghosts in the scalar sector at high momenta, and are thus quantum-mechanically unstable 

m- 

Let us conclude by emphasizing that the stability conditions we have derived in this 
Appendix are necessary but not sufficient. Because we have not analyzed the stability of all 
Fourier modes, a violation of our conditions points to the instability of certain field modes, 
but their satisfaction does not imply that all the modes are free of instabilities. 
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